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Abstract 


In this paper, we consider a special case of Weber location problem which 
we call goal location problem. The Weber location problem asks to find 
location of a point in the plane such that the sum of weighted distances 
between this point and n existing points is minimized. In the goal location 
problem each existing point P; has a relevant radius r; and it’s ideal for us to 
locate a new facility on the distance r; from P; for i = 1,...,n. Since in the 
most instances there does not exist the location of a new facility such that its 
distance to each point P; be exactly equal to r;. So we try to minimize the 
sum of the weighted square errors. We consider the case that the distances 
in the plane are measured by the Euclidean norm. We propose a Weiszfeld 
like algorithm for solving the problem and also we use two modifications of 
particle swarm optimization method for solving this problem. Finally the 
results of these algorithms are compared with results of BSSS algorithm. 


Keywords: Location theory; Weiszfeld method; Particle swarm optimiza- 
tion. 


1 Introduction 


Location theory is onc of the uscful and interesting ficlds of operations re- 
search and has many applications in real life. Single facility location problems 
are fundamental problems in location theory in which n existing demand 
points are located at known distinct points P,,--- ,P,, and a new facility 
should be located at a point X such that the sum of weighted distances from 
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it to demand points is minimized. Let d(X, P;) represent the distance trav- 
eled per trip between points X and P;, the weight w; represents the product 
of cost per unit distance traveled and number of trips made per year between 
the new facility and existing demand point i and some times shows value of 
demand in this point. The total cost is given by: 


min F(X) = ys wd(X, P,). (1) 


The Euclidean problem is obviously referred to as the Steiner-Weber problem 
or the general Fermat problem, and has an extraordinary longevity (see [8]). 
In fact, a version of the problem for the case n = 3 and w; = 1 for i = 1, 2,8, 
was posed purely as a problem in geometry, by Fermat early in the seven- 
teenth century, and was solved by Toricelli prior to 1640 (see [16]). The 
problem was studied by Steiner, a Swiss mathematician, in the nineteenth 
century, and by Weber, a German economist, early in the twentieth cen- 
tury. Single facility location with Euclidean distance is also called Weber 
problem and an iterative procedure was proposed for solving this problem by 
Weiszfeld [24] in 1937 and rediscovered by Kuhn [17] in 1962. This algorithm 
is widely used because of its simplicity and effectiveness. The algorithm can 
be generalized to other location problems where the cost is a function of the 
Euclidean distance rather than just being proportional to the distance (see 
e.g. [10]). Some other cases of Weber location problem can be find in [4,5,23]). 


Fathali et al. [12], considered a special case of Weber problem in which 
a specified radius for every demand point is considered and the distance be- 
tween the new facility and the demand point P; is equal to the corresponding 
radius, r;. However since in the real instances rarely exist the location of a 
new facility such that its distance to each point p; is exactly r;. So they tried 
to minimize the sum of the weighted square errors. We call this problem as 
Goal Square Weber Location Problem (GSWLP). Jamalian and Fathali [14] 
presented a linear model for the location problem with minimum absolute 
error on the block norms. 


In what follows we explain the GSWLP in Section 2. Model properties 
and results are stated in Section 3. In this section main results and a proposed 
algorithm for solving problem are given. Because the objective function may 
have many local optima, we also use particle swarm optimization method for 
solving this problem. In Section 4, the particle swarm optimization method 
is explained. In Section 5, the computational results of comparing four algo- 
rithms WLA, PSO, PSOC and BSSS for this problem are given. 
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2 Goal Square Weber location problem 


In single facility location problem, we want to find location of a new facility 
such that the sum of weighted distances from the facility to all demand 
points is minimized. In Fathali et al. [12] a radius, r;, and value of demand, 
w; correspond to every point P; is considered. We want to find the location 
of a new facility such that the distances between it and the demand points 
P;, i=1,...,n is equal to r;. Since this situation occurs in real applications 
rarely we try to minimize sum of weighted squared error. 

In fact, in this problem for every point a goal distance is considered and we 
want to minimize the error of satisfaction location of facility in these goals. It 
is desirable for us that new facility be located on circle of correspond radius 
and center of demand point. For an example of this problem, if r; = r, 
i = 1,2,---,n and demand points be located on a circle of radius r, the 
optimal point is center of circle. In the case r; = r, i = 1,2,---,n the 
problem converted to the square Weber problem, so we call this problem as 
Goal Square Weber Problem. 

The problem studied in this paper is different from the covering location 
problem (CLP). Two important covering models are the maximal and min- 
imal covering location problems which are suitable for siting desirable and 
undesirable facilities, respectively (see e.g., [8,6]). In both of these problems a 
facility covers a demand point if the demand point lies within a pre-specified 
coverage radius. The objective is to locate a number of facilities so as to min- 
imize or maximize the total coverage. There are some differences between 
CLP and GSWLP. Firstly, CLP is a multifacility location problem, while in 
the GSWLP we want to find the location of a single facility. Secondly, in 
CLP a radius is associated with facilities, in contrast to the GSWLP, where 
a radius r; is associated with demand point P;. Thirdly, in the maximal 
CLP, it is desirable if the distance between the demand points and facilities 
be less than or equal to a pre-determined radius and in the minimal CLP this 
distance should be more than or equal to a given radius, while in GSWLP it 
is desirable if the new facility exactly be in the given distances of the demand 
points. 

Determination of the location of a company in the vicinities of some cities 
such that the setting up costs and the transportation costs and minimized can 
be an application of this problem as mentioned in [12] . Suppose that the cost 
of establishing a facility in the regions that are farther than a given distance 
r,; from city r; is very low. On the other hand a move away from a city causes 
the transportation costs to increase. Therefore a tradeoff between the setting 
up costs and the transportation costs seems to be reasonable. GSWLP also 
has some other applications for locating facilities that are either desirable or 
undesirable. In these cases we want the facility not be close than a specified 
distance to facility centers, because of its undesirability, on the other hand, 
if the facility be so far from the facility centers, cost of providing security, 
human forces, transportation installation, and other costs will increase. 
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3 Model Properties 


Let P; = (a;,0;) for i = 1,--- ,n be given points in the plane. Formulation 
of GSWLP is as follows: 


min F(X) = » wi (d(X, P;) — ri)? (2) 


where X = (2,y), d(X,P;) = \/(a—a;)? + (y—6;)?, and rj and w; are 
corresponding radius and weight of demand point P;, respectively. Fathali 
et al. [12] show that the objective function of model (2), is non convex and 
optimal solution is in extended rectangular hull of demand points. Extended 
rectangular hull of demand points is the smallest rectangle which contains 
all the demand points with their circles. We present a new proof by Farkas 
lemma in the Theorem 1. The advantage for our new proof is that in this 
proof we show there is an improvement direction from any point which is out 
of rectangular hull. 

Farkas lemma is one of the theorems of alternatives a good discussion of 
which is given by Mangasarian [19]. It can also be found in Dantzig and 
Thapa [9] and Bazaraa et al. [2], among others. 


Lemma 1. [Farkas Lemma (1902)] For a given matric A and a vec- 
tor C, either AX > 0 and CX < 0 has a solution, or there exists a vector 
like W such thatWA=C andW > 0. 


Theorem 1. The optimal solution of Equation (2) is in the extended rect- 
angular hull of demand points. 


Proof. First we define four points as follow: 


RHA, = min {a; — r;|i = 1,2,--- ,n} 
RH» = max {bj + r;|i = 1,2,--- ,n} 
RH3 = max {a; + rj|t = 1,2,--- ,n} 
RHA, = min {b; — r;|i = 1,2,--- ,n}. 


We show that the optimal point is in the rectangular hull of these four points. 
By contradiction, we suppose that X, = (ap,yp)" is the optimal point of 
model (2) and is out of rectangle. Therefor it is out of cone generated either by 
(RH, — RH3,0) and (0, RH4— RHA2) or (RH3— RH;,0) and (0, RH — RH4). 
Case 1: suppose that the optimal point is out of cone generated by (RH, — 
RH3,0) and (0, RH4— RH2). So the vector X, —(RH3, RH2)" is not in the 
cone, i.e. following system has not solution: 
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vi (RA, = RHs3,0) + v9 (0, RHA, - RH) = (xp — RA3, yo — RH2) 
V1, V2 = 0. 


By Farkas’ lemma, the following system have solution: 


RHA, — RH3 0 U1 
0 RHA, ca RHA» U2 


[xp _ RA3, yp _— RH») Ei <0. 
ug 
So by the first inequality of (3) we have: 


(RA; foe RH3)uy > 0 
(RH, = RH2)u2 > 0 


since RH, < RH3 and RH, < RH» therefore uj < 0 and ug < 0. And by 
the second inequality of (3) we obtain 


(xp _ RH3)uy + (Yo = RH2)u2 <0. 
Moreover a; < RH3 and b; < RH» for i = 1,2,--- ,n so 
(xp _ aj )Uy + (Yo _ bj )ue < (xp _ RH3)u1 + (Yo = RH2)u2 <0. 
Now, if z» — RH3 < 0 and y— RH2 > 0 then let U = (2) and if 
2 


U1 


ty ~ REly > O and yy RHfp <0 then tet 0 = (“4 


) ana if 25 — Rey > 0 


and y, — RH2 > 0 then let U = eS . So in each case 
2 


Therefore 
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This means that U is a descent direction and X» Con not be optimal point. 


pote that in the last equality statement 1 — UX Py is positive for all 
y because X, is out of extended rectangular hull and d(Xy,P;) > 1; so 
UX PD py Sl. 
oP) 


Case 2: suppose that optimal point is out of cone generated by (RH3 — 
RH;,,0) and (0, RH. — RH,). Proof of this case is the same as Case 1. So 
the optimal point is in rectangular hull generated by RH), RH2, RH3 and 
RAHy4. 


4 Weiszfeld-like algorithm 


In this section, we propose an iterative algorithm for solving GSWLP. Our 
algorithm is based on the Weiszfeld method. The original Weiszfeld algorithm 
is an iterative method for solving the following Weber problem. 


Mink (X = Yom (2 — ai)? + (y — 6;)?. (4) 


aF(X*) 


From the necessary condition of optimality = 0 for this problem, 
the Weiszfeld method starts with an initial solution and try to find optimal 
solution by the following iterative relations; 


sit 


vis 1 Vices ai)? +(yn—b;)? 


Lk+1 n 
Dist Toran? torch? 


n w;7o; 
= Diml Toray? 
Yk+1 > 3H Wi 7 


1 (ax —ai)? +(yx—bi)? 


Now consider the GSWLP. The same as Weiszfeld method from the nec- 
essary optimality condition to obtain a candidate for optimal solution we 


set 
ae = 2d (1 ER) P,) =0 (5) 


therefore 

Dia Wi (1 7 ner) f 
ie Wi (1 = er] 

Like in the Weiszfeld algorithm, there is a possibility that during the 


iterations, d(X, P;) is zero for one of the points, and then the ratio is infinite 
and the process cannot be continued. For practical purposes, the whole issue 


X= (6) 
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of singular points may be avoided by replacing each distance in (6) by its 
hyperbolic approximation (e.g., see [18]) given by: 


Tie (1 — ayes) Pi 
naw (1- aretha) 


where € > 0 is a smoothing constant. Therefor we can use the following 
iterative method. 


X= (7) 


Algorithm [WLA]. 
Initialization: 


1. Choose an initial solution, Xo and set k = 0. 


Iteration step: 
While not stopping criterion do 


1. 
Dini Wi (1 Sse) = 
Tk+1 = 
ier Wi (1 ee 
Viet Wi (1 ars) i 
Yk+1 = 
int Wi (1 contig ay) 
2. k:=k+1 
endwhile 


The stopping criterion rule is usually one of the following criteria: maxi- 
mum number of iterations or reaching a tolerance. 

Now we discuss the convergency of WLA. From Equation (6) if we con- 
sider the right hand side of this equation as X, and the left hand side of it 
as X,41, So we have 


ae Wi (1 = ao) P; 
Diat Vi (1 = avers) 


then by choosing an initial solution, Xo, for k = 1,2,--- we have 


(8) 


Xey1 = 


72 J. Fathali and A. Jamalian 


1 z Ti; 
X41 = Xk Sow (1 % ) om P;) 
in Wi (1 Z a) i=1 (Xx, Pi) 
(9) 
Let us assume that 
Hesse > (41 = (=P) (10) 
Ke QZ Wi d(X;,, P;) k a 


it is clear that dz is a descent direction, i.e. dy VF(Xx) < 0, so the Equation 
(9) is a gradient method and the sequence {F'(X;,)} is non increasing. Let 


1 OF (X) 
Im Wi (1 = ans) os 


The function F(X) is two times differentiable and h(X) in X 4 P,, 
i= 1,--- ,n, is continuous. 


h(X) =X (11) 


Lemma 2. The sequence generated by Equation (8) is in affine hull of the 
Pit = 1,2,--- Nn. 
Doin Pi 


Proof. In each iteration, we have Xx41 = rE? 
c= Te 


, in which 
Oo; = ath (1 - aep5): Let % = oak aj we have 7; € [-1,1] and 
1% = 1. From Equation (8) the sequence {X;} is a affine combination of 
demand points. So the sequence {X;,} is in affine hull of P;,2 = 1,2,---n. 


Lemma 3. If X_ = Xx41 then OFS) =0. 


Proof. From the Equation (9), we have Xx41 = h(X;,) = Xp - By ees) 
where Ay = 1____. If Xnii = Xp then —2d, 228 = 0, ie. 
Thi wi (1- aes) 
dF(Xx) _ g 
OX — ™ 


Lemma 4. If X;, x Xp then F(Xp41) < F(X), and If F(Xx) = P(X ) 
then Xk — Xk41- 


Proof. From the Equation (9), if X;, 4 X41, we have OF Se) #0. Let 
D = —2EQ) 59 DT.VF(Xx) < Oie. F(Xn41) < F(Xz). 

To show the second part of lemma if F'(X;,) = F(X;,41) then since for every 
k, the direction dy, is non increasing we have X, = Xx41. 


Note that the direction d, in iteration k, is a descent direction, so the 
sequence {F'(X;,)} is non increasing. Since this sequence is bounded from 
below by zero, the limit limp... F(X,) exists. 
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Theorem 2. The sequence {X,} converges to a local minimum or a saddle 
point of function F(X). 


Proof. By Lemma 2, the sequence {X;,} is in affine hull of demand points 
and since by explanation of Equation (9), df are descent direction so the 
sequence {X;,} is converged to a point like X. Since limg_,.. F(X;) exists, 
limpsoo F(Xk41) — F(Xz) = 0 and F(h(X)) = F(X). By Lemma 4, we 
have h(X) =X. If Xp41 # Xz then the inequality F(X,41) < F(X) holds 
strictly. From continuity of function h(X), if X, tends to X then X,41 
will tend to h(X) = X and {X;} converges to X. According to Lemma 3, 


gues = 0 and since {F'(X;,)} is non increasing and bounded from below 


by zero, the point X con not be local maxima, in other word, it is a local 
minimum or saddle point of function F(X). 


We presented an analytic method for solving GSWLP that do not guar- 
antee converging to a minimum and may fall in saddle point or local minima 
which is no minimum. Since the objective function of GSWLP is non convex, 
the local minimum may not be global (see Example 1). So we should use a 
global optimization method. In the next section we apply a particle swarm 
optimization algorithm to find the best solution of problem. 


Example 1. Consider the points, their weights and their radiuses that are 
given in Table 1. This problem is solved in [12] by big square small square 
method and the value of optimal solution is 182. But if we solve this problem 
by WLA we will see this method convergence to a point which its value of 
objective function is 209. 


Table 1: Data for 18 points problem. 


b 


2 


WwWNwrnwwls 
FPmnNmwWN bd 


5} 
3 
5 
6 
8 
1 


Note that this is an special example that WLA traps in a local minima. 
This case rarely happens. We examined many other test problems but we 
could not find any example with this property. Also as we shall see in the 
computational experiments section, in non of the examined test problems in 
this section this case dose not happen. 
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5 Particle Swarm Optimization 


Nature has always been an inspiration for researchers in designing solutions 
for many problems. Particle Swarm Optimization (PSO) is a meta-heuristic 
algorithm that is inspired of the behaviors of social models like bird flocking 
or fish schooling. It is introduced by Eberhart and Kennedy [11], as an 
optimization method for continuous nonlinear functions. Later, it has been 
applied to wide range of problems due to its conceptual and implementation 
simplicity (see e.g. [1,15]). 


Particle swarm optimization is a stochastic optimization technique based 
on individual improvement, social cooperation and competition in a popula- 
tion. PSO is famous for its fast convergence to a solution close to the optimal, 
by balancing the local search and global search i.e., exploitation and explo- 
ration, respectively. PSO is a population based method in which a swarm 
includes several individuals called particles. Each particle has a position and 
velocity vector. Particles are initially positioned randomly. The position of 
a particle is a candidate solution and it is updated at each iteration by using 
the particle’s current velocity. The velocity of a particle is updated by using 
the particle’s inertia weight as well as the social interaction (its tendency to 
follow the best direction experienced by the group) and personal experience 
(the best direction experienced by itself) of the particle. In the course of 
several iterations, particles make use of this experience and are supposed to 
move towards the optimum position. 


The stopping criterion rule is usually one of the following criteria: maxi- 
mum number of iterations, the maximum CPU time, the number of successive 
best objective function values without any improvements or reaching a pre- 
determined tolerance. During the update step, for each particle, the velocity 
and position vector of the particle at iteration t+ 1 are calculated by using 
the Equations (12) and (13). 


ai — Wr} 5 +11 (pbest; , _ v5) + cere (gbest’ — x; ) (12) 
cs = xs + Wie (13) 


t 
inj 
dimension of the i” particle at iteration t, respectively. The parameters c, 


and c2 are coefficients of learning factors, which are the weights of contribu- 
tions of personal experience and social interaction. The stochastic behavior 
of PSO is achieved by 71 and re which are random numbers, generally in 
(0,1). The parameter w is the inertia weight which is the balancing factor 
between exploration and exploitation. Small inertia weight facilitates more 
exploitation while large inertia weight enables more exploration. 


In these equations, vj, and x}, are the velocity and position of the gee 


After the update step, the fitness function value is calculated for each par- 
ticle based on its position. The local best position pbest of each particle and 
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the global best position gbest are updated using these fitness values as follow: 


t41 _ f pbestt F(att') > F(pbest’) 
ae Olle { ait} F(zit1) > F(pbest?) 14) 
gbest’** = argminies,2,... ny F (pest; **). (15) 


These ideas lead to the following algorithm. 
Algorithm [PSO]. 


Initialization: 
1. Initialize iteration counter. 
2. Initialize N random position of particles and store them in S. 
3. Initialize N random velocities and store them in V. 
4. Initialize N pbest and store them in P. 
5. Set gbest’ equal the best pbest in P. 


Iteration step: 
While not stopping criterion do 


1. For each i*” particle: 


a) Update V: Calculate velocity v'*! using (12 
y 4: gs 


(b) Update S: Calculate position «/*? using (13) 
t41 


4 


(c) Update P: Calculate position pbest;"” using (14) 


2. Update gbest: gbest't? = arg minjes19.... Ny { f (pbesti**)} 
3. t:=t+1 


endwhile 


Where V is the set of velocities and P is the set of pbests. To apply this 
algorithm for Goal Square Weber Location Problem we set S' the region of 
optimal solution, i.e. extended rectangular hull of points and each particle 
as a solution. In this algorithm first we choose N random positions and 
velocities for each particle. From rules of particle swarm optimization we 
obtain the pbest and gbest of swarm. The algorithm with updating variables 
continues until termination conditions reached. 


Many modification for PSO algorithm has been proposed to improve the 
performance of this algorithm. In standard PSO algorithm, the information 
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of individual best and global best are shared by next generation particles. 
Shi and Eberhart [21] use dynamic inertia weight that decreases according 
to iterative generation increasing. The weight w in (12) follows the relation 
Wit1 = aw, where 0 < a < 1. A large inertia weight facilitates global 
exploration and a smaller inertia weight tends to facilitate local exploration 
to search the local neighborhoods in the current search area more effectively. 

Clerc and Kennedy [7] present another good modification of PSO. They 
consider the form of a constriction coefficient kK which controls all the three 
components in velocity update rule. This has an effect of reducing the velocity 
as the search progress. In this modification, the velocity update is given as: 


ee =k (vi. 5 + e171 (pbest} , = xis) + core (gbest’ = xi ;)) (16) 


where 


2 
K= 17 
|2-¢-— /¢? — 4 | a 
where @ = cy + cg > 4. This version is referred to as PSO with constriction 
or PSOC. We apply the original PSO and PSOC for solving GSWLP. 

In the next section we compare the WLA with PSO, PSO with constric- 
tion (PSOC) and big square small square (BSSS) method. Big Square Small 
Square (BSSS) method is a geometrical branch and bound algorithm origi- 
nally suggested by Hansen et al. [13] for solving an obnoxious facility location 
problem. The idea is to divide the plane into regions (squares), over each 
of which a lower (upper) bound for the problem is found. If the objective 
over a given square is worse than an existing upper (lower) bound, then that 
square is fathomed. The procedure continues until a prescribed tolerance is 
achieved. This method also applied by some authors such as Plastria [20], 
Skriver and Andersen [22], Zaferanieh et al. [25] and Fathali et al. [12] for 
solving location problems. 


6 Computational results 


The proposed algorithms were tested on forty test problems with 100 to 
3000 points that are generated randomly in MATLAB. In these problems 
demand point’s coordinates, weights and radiuses, all are positive numbers 
and generated randomly in [1,60], [1,3] and [1, 10], respectively. 

We compare the Weiszfeld like algorithm (WLA), PSOC and PSO with 
BSSS method of Fathali et al. [12]. Algorithms are coded in MATLAB soft- 
ware and run on a Laptop with Core(TM)2 CPU with 2.00 GHz processor 
and 1 GB of RAM. In the PSO and PSOC algorithms we should generate 
randomly N initial position of particles and velocities. Therefore in differ- 
ent executions of these methods we may find different solutions, so we run 5 
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times this two methods for all problems and report the average results. In all 
methods the iteration step of the algorithm is repeated until the tolerance is 
achieved. The tolerance for all examples was taken to be 0.01. 

In our experiments with different values of cy and c2 (c1,c2 € [1,3]) we 
found the best values to be cy = cg = 2.1. The inertia weight is started 
with w = 1.5 and decrease in each iteration such that reach 0.2 at the end 
of algorithm. The population size is set N = 50 for problems with n = 
100,--- ,500 and N = 100 for remaining problems. We also examined larger 
numbers of population size. Our tests showed that larger population size 
slow down the algorithm but did not lead to better solutions. 

Table 2 shows the computational results. The column with the heading 
%Error indicates the relative error; i.e., 


f = Soest 


| foest 


| - 1000 


where foest is the best solution obtained by the four methods and f is the 
value of objective function obtained by each method. 

Note that the BSSS method continues until the side of the sub-square is 
less than a given tolerance. So the solutions which obtained by this method is 
near optimal with the given tolerance. By this method we can also compute 
a lower bound for the objective for each sub-square. The column with the 
heading LB in Table 2 indicates the lower bound for the objective for the 
best sub-square. This sub-square contains optimal solution. Also the column 
with the heading gap shows the relative gap between the lower bound and 
the best solution, i.e., 

) Soest —LB 


- 1000. 
foest | 


We observe the largest relative gap is 5.79 x 107+, which happens for the test 
problem number 23. 

Table 3 contains the average total CPU times for test problems. The 
CPU times of WLA are less than other methods in all cases. However the 
other methods could find better solution in some test problems. 


To compare speed of presented methods one may consider the time com- 
plexity of these methods. However since these methods are improvement 
methods then with more iterations we may obtain better solutions. There- 
fore in the following we give their time complexity per iteration. 

The main step of WLA is computing x and y in each iteration. Since we 
need 23n operations for each of x and y, obviously, the time complexity per 
iteration is O(n). Where n is number of existing points. 

In step 1 of each iteration of PSO, we should calculate velocity v, position 
of « and pbest according to (12), (13) and (14), respectively. They need 
O(nN) operations. Where N is the population size. Step 2 of Iteration step 
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Table 2: The objective functions and errors 
objectivefunction Error 
NO. n LB WLA BSSS PSO PSOG WLA BSSS PSO _PSOG gap 
T 100 59599 59617 59618 59617 59617 0.000 0.017. 0.000 0.000 | 0.302 
2 100 68958 68995 68995 68995 68995 0.000 0.000 0.000 0.000 | 0.536 
3 100 82265 82287 82287 82288 82287 0.000 0.000 0.012 0.000 | 0.267 
4 100 78156 78232 78180 78182 78180 0.665 0.000 0.026 0.000 | 0.435 
5 100 76104 76145 76145 76145 76145 0.000 0.000 0.000 0.000 | 0.539 
6 200 121305 121340 121341 121341 121340 | 0.000 0.008 0.008 0.000 | 0.288 
7 200 166160 166254 166254 166255 166254 | 0.000 0.000 0.006 0.000 | 0.565 
8 200 159632 159720 159720 159720 159721 | 0.000 0.000 0.000 0.006 | 0.551 
9 200 162253 162312 162312 162312 162312 | 0.000 0.000 0.000 0.000 | 0.363 
10 200 160276 160324 160324 160325 160324 | 0.000 0.000 0.006 0.000 | 0.299 
11 300 233660 233765 233766 233765 233766 | 0.000 0.004 0.000 0.004 | 0.449 
12 300 234183 234314 234314 234314 234315 | 0.000 0.000 0.000 0.004 | 0.559 
13 300 245908 246045 246045 246046 246045 | 0.000 0.000 0.004 0.000 | 0.557 
14 300 241698 241761 241761 241761 241761 | 0.000 0.000 0.000 0.000 | 0.267 
15 300 235575 235674 235621 235642 235621 | 0.225 0.000 0.089 0.000 | 0.195 
16 400 309670 309777 309777 309777 309776 | 0.003 0.003 0.003 0.000 | 0.342 
17 400 337358 337543 337544 337543 337543 | 0.000 0.003 0.000 0.000 | 0.548 
18 400 319847 320030 320030 320031 320032 | 0.000 0.000 0.003 0.003 | 0.493 
19 400 326612 326734 326735 326734 326735 | 0.000 0.003 0.000 0.003 | 0.373 
20 400 318651 318783 318783 318783 318783 | 0.000 0.000 0.000 0.000 | 0.389 
21 500 394201 394364 394364 394364 394364 | 0.000 0.000 0.000 0.000 | 0.418 
22 500 391482 391654 391598 391597 391598 | 0.143 0.003 0.000 0.003 | 0.294 
23 500 407515 407751 407752 407751 407751 | 0.000 0.003 0.000 0.000 | 0.579 
24 500 412614 412763 412763 412763 412763 | 0.000 0.000 0.000 0.000 | 0.361 
25 500 402198 402348 402277 402277 402277 | 0.176 0.000 0.000 0.000 | 0.196 
26 1000 | 786520 786731 786732 786732 786732 | 0.000 0.001 0.001 0.001 | 0.268 
27 1000 | 796696 797148 797149 797148 797148 | 0.000 0.001 0.000 0.000 | 0.191 
28 1000 | 812939 813407 813406 813407 813406 | 0.001 0.000 0.001 0.000 | 0.575 
29 1000 | 805501 805634 805634 805634 805634 | 0.000 0.000 0.000 0.000 | 0.165 
30 1000 | 802280 802431 802432 802431 802432 | 0.000 0.001 0.000 0.001 | 0.313 
31 2000 | 1465017 1465731 1465705 +=—-:1465718 + «1465705 | 0.018 0.000 0.009 0.000 | 0.469 
32 2000 | 1582076 1582976 1582977 1582977 1582977 | 0.000 0.001 0.001 0.001 | 0.569 
33 2000 | 1633068 1634000 1634000 1634000 1634000 | 0.000 0.000 0.000 0.000 | 0.570 
34 2000 | 1563967 1564236 1564185 1564186 1564185 | 0.001 0.000 0.001 0.000 | 0.139 
35 2000 | 1607103 1607640 1607640 1607641 1607640 | 0.000 0.000 0.001 0.000 | 0.334 
36 3000 | 2297871 2298872 2298873 2298872 2298873 | 0.000 0.001 0.000 0.001 | 0.435 
37 3000 | 2362318 2363663 2363663 2363663 2363663 | 0.000 0.000 0.000 0.000 | 0.569 
38 3000 | 2403658 2405041 2405040 2405040 2405040 | 0.001 0.000 0.000 0.000 | 0.575 
39 3000 | 2387987 2388745 2388746 2388745 2388753 | 0.000 0.001 0.000 0.003 | 0.317 
40 3000 | 2324214 2325395 —«-2325396 «2325396 —«2325395 | 0.000 0.001 _0.001_~_—i0.000_-| _0.508 
Total 1.233 0.051 0.172 0.030 


Table 3: The CPU Time of PSOC, PSO, WLA and BSSS algorithms 


n WLA ~ BSSS PSO = PSOC 
100 | 0.002 4.204 0.630 0.778 
200 | 0.004 7.765 1.145 1.318 
300 | 0.005 11.363 1.686 1.852 
400 | 0.007 15.061 2.192 2.920 
500 | 0.008 18.205 2.709 3.470 
1000 | 0.014 83.607 5.348 5.929 
2000 | 0.027 389.972 10.570 12.346 
3000 | 0.040 576.231 15.750 17.230 


Efficient methods for goal square Weber location problem 79 


Table 4: The objective functions and errors in a constant time 


objective function %Error 

test# n time BSSS PSO PSOCG BSSS PSO PSOCG 
T 100 0.5 59630 59617 59617 0.218 0.000 0.000 
2 100 0.5 68997 68996 68995 0.028 0.014 0.000 
3 100 0.5 82290 82289 82287 0.036 0.024 0.000 
4 100 0.5 78241 78205 78202 0.780 0.320 0.281 
5 100 0.5 76147 76146 76145 0.026 0.013 0.000 
6 200 1.0 121355 121345 121343 0.124 0.041 0.025 
7 200 1.0 166799 166263 166260 3.278 0.054 0.036 
8 200 1.0 160699 159771 159730 6.129 0.319 0.063 
9 200 1.0 162330 162318 162314 0.111 0.037 0.012 
10 200 1.0 160345 160325 160324 0.131 0.006 0.000 
11 300 1.5 233809 233766 233767 0.188 0.004 0.008 
12 300 1.5 234316 234316 234316 0.008 0.008 0.008 
13 300 1.5 246072 246047 246045 0.110 0.008 0.000 
14 300 1.5 241763 241770 241761 0.008 0.037 0.000 
15 300 1.5 235704 235664 235651 0.352 0.182 0.127 
16 400 2.0 309790 309779 309780 0.039 0.009 0.012 
17 400 2.0 337568 337543 337543 0.075 0.000 0.000 
18 400 2.0 320116 320090 320076 0.258 0.180 0.138 
19 400 2.0 326814 326784 326784 0.248 0.153 0.153 
20 400 2.0 318853 318783 318783 0.220 0.000 0.000 
21 500 2.5 394384 394366 394365 0.061 0.005 0.003 
22 500 2.5 391755 391598 391599 0.403 0.003 0.005 
23 500 2.5 407801 407751 407751 0.151 0.000 0.000 
24 500 2.5 412833 412772 412763 0.211 0.027 0.000 
25 500 2.5 402451 402279 402280 0.433 0.006 0.009 
26 1000 5.0 786851 786742 786733 0.153 0.011 0.002 
27 1000 5.0 797153 797149 797149 0.004 0.001 0.001 
28 1000 5.0 813407 813410 813407 0.001 0.004 0.001 
29 1000 5.0 805710 805642 805639 0.094 0.008 0.005 
30 1000 5.0 802671 802431 802433 0.299 0.000 0.002 
31 2000 10.0 1465978 1465720 = 1465711 0.186 0.010 0.004 
32 2000 10.0 1583288 1582977 1582977 0.197 0.001 0.001 
33 2000 10.0 1634121 1634004 1634000 0.121 0.004 0.000 
34 2000 10.0 1564480 1564186 1564186 0.189 0.001 0.001 
35 2000 10.0 1607871 1607665 1607660 0.144 0.025 0.020 
36 3000 15.0 2299001 2298876 2298874 0.056 0.004 0.002 
37 3000 15.0 2363820 2363663 2363663 0.066 0.000 0.000 
38 3000 15.0 2405131 2405041 2405041 0.037 0.001 0.000 
39 3000 15.0 2388876 2388754 2388759 0.055 0.004 0.006 
40 3000 15.0 2325480 2325396 —«2325395 0.037 0.001 0.000 
Total 15.265 1.525 T.04 


calculate gbest which needs O(V) operations. Therefor per iteration of PSO 
has O(nN) time complexity. 

PSOC is the same as PSO just we need calculate velocity v by (16). It also 
need O(nN) operations. Therefor the time complexity of PSOC is O(nN) 
for per iteration. 

Due to [12] in each iteration of BSSS we should calculate the objective 
function for a given point in a square. This take O(n”) time, therefor each 
iteration of BSSS has O(n”) time complexity. 

To compare the algorithms on equal terms, we allow the BSSS, PSO and 
PSOC algorithms run during the same time. Table 4 contains computational 
results for this case. Since WLA is very fast, its results for the given times 
in Table 4 is the same as Table 2. 


7 Summary and conclusion 


In this paper we studied a new version of the single facility location problem in 
which we want to find location of a new facility such that the sum of weighted 
squared errors is minimized. In general, this problem is non convex and we 
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showed that optimal solution of problem is in extended rectangular hull of 
demand points. We proposed two approaches for solving this problem. An 
iterative Weiszfeld-like procedure and two modifications of PSO algorithm. 
We compared the results with those obtained by BSSS algorithm. It was 
shown that for almost all problems the PSOC outperforms the other three 
approaches. 
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